Emergence of Wigner molecules in one-dimensional systems of repulsive fermions 

under harmonic confinement 
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A Bethe-Ansatz spin-density functional approach is developed to evaluate the ground-state den- 
sity profile in a system of repulsively interacting spin- 1/2 fermions inside a quasi-one-dimensional 
harmonic well. The approach allows for the formation of antiferromagnetic quasi-order with increas- 
ing coupling strength and reproduces with high accuracy the exact solution that is available for the 
two-fermion system. 
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Introduction — It is known from bosonization [l| that 
the ground-state |GS) of a one-dimensional {ID) sys- 
tem of repulsively interacting spin-1/2 fermions carries 
antiferromagnetic quasi-ordering. Quantum fluctuations 
prevent any type of true spin-symmetry breaking in 1_D 
and |GS) may be viewed as a superposition of two spin- 
density waves (SDW), 



iGs)cxinn-) + UTiT-) 
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so that the local spin polarization is rigorously zero. Nev- 
ertheless, ordering is signaled by the co-called ^^2kp 
Akp crossover" in the wave number q of bulk Friedel oscil- 
lations: for vanishing and weak interactions the ground 
state is liquid-like with q = 2fcF, whereas at strong re- 
pulsions a periodicity q = dfcp emerges. Here fcp is the 
Fermi wave number, related to the fermion density n by 
fcp — 7rn/2, and each modulation is weighted by a slow 
decay factor. 

In previous work Q we studied a system of N spin-1/2 
fermions subject to a strongly anisotropic harmonic po- 
tential, characterized by angular frequencies uj^_ and cij|| 
in the radial and axial directions with wy <C cjj^. The 
fermions are dynamically ID if the anisotropy parame- 
ter of the confinement is much smaller than the inverse 
particle number, o^n/o;^ ^ N^^. The two species of 
fermionic particles are taken to have the same mass m, 
to be in equal numbers N-^ and , and to interact via a 
contact repulsion with an effective ID coupling strength 
5iD > [3| . The system is described by the Hamiltonian 



unit of length and the harmonic-oscillator quantum hiu^^ 
as unit of energy, the Hamiltonian ^ can be shown to 
be governed by the dimensionless coupling parameter 
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For Vcxt — the Hamiltonian ^ reduces to the homoge- 
neous Gaudin-Yang model, which is exactly solvable by 
means of the Jiethe-Ansatz technique 4]. 

In Ref . i2 ^ density- functional scheme using as reference 
the Bethe- Ansatz homogeneous fluid was proposed to cal- 
culate the ground-state properties of the inhomogeneous 
ID system described by the Hamiltonian However, 
this approach was found to fail at strong-coupling (see 
Fig. 8 at A = 10 in Ref.i3): more precisely, the exchange- 
correlation potential proposed there is unable to describe 
the 2fcF dfcp crossover, which as noted above is ex- 
pected to occur upon increasing the strength of the re- 
pulsive interactions between antiparallel-spin particles. 
In the present note we propose a simple functional that 
embodies this crossover and is capable of describing in- 
homogeneous Luttinger liquids with strong repulsions. 
The main idea is to capture the tendency to antifer- 
romagnetic spin ordering (i) by adding an infinitesimal 
spin-symmetry-breaking magnetic field to the Hamilto- 
nian (0), 
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where contact interactions between parallel-spin parti- 
cles are suppressed by the Pauli principle and Vcxt = 



J2i=i^oxt{zi) = {'muj'^^/2)J2i=i zf is the external static 
potential associated with the axial confinement. Choos- 
ing the harmonic-oscillator length a|| = y^h/{mLu^^) as 
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with iiB the Bohr magneton; and (ii) by resorting to 
spin- density functional theory (SDFT) 

The ground-state energy of the Gaudin- Yang model 
for finite {Sz) — The addition of a magnetic field to 
the Hamiltonian ^ as prescribed by Eq. ^ requires 
that we know the ground-state energy of the homoge- 
neous Gaudin-Yang model for the situation ^ N^, 
which will be our reference system with Luttinger-liquid 
ground-state correlations. In the thermodynamic limit 
{N-^jNi, and L ^ oo, L being the system size) the 
properties of the reference fluid are determined by the 
linear density n = N/L, by the spin polarization ( = 
(iV| — Ni)/N, and by the effective coupling giD. The 



2 



linear density and the coupling parameter can be conve- 
niently combined into a single dimcnsionless parameter 
7 = rngin/ifi^n). 

The energy per particle can be written in terms of the 
momentum distribution p{k) as 
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The function p(k) can be calculated by solving the 
two-coupled Gaudin-Yang Bethe-^nsate integral equa- 
tions (ll, 
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and 
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where Q and _B are determined by the normalization con- 
ditions 



p{k)dk = 



and 



f\ik)dk = ^{1-0. 
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We have fitted the results for ecsin, Ci 9id), obtained 
from the numerical solution of the Gaudin-Yang equa- 
tions ([6])- ([9]), with a simple yet accurate parametrization 
formula given below. 

We introduce the Fermi wave number kp = 7m/2 of the 
unpolarized system and the corresponding Fermi energy 
£p — h'^kp/{2m). We write 



£Gs{n, C, 5id) = K{n, C) + Secsin, C, 5id) 



where 
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is the kinetic energy of the noninteracting system per par- 
ticle. We find that &gs(^:Cj5id) in units of the Fermi 
energy, / = fcGs(fi, d 5'iD)/eF, can be very accurately 
represented by the parametrization 
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[e{x) ~l/3]{l + aix)C^ 
- [l + a(x)+/3(x)]C^} (12) 
where x = e{x) is given in Eq. (21) of Ref.'?, and 
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FIG. 1: The interaction contribution / to the ground-state 
energy Secsin, giu) of the homogeneous Gaudin-Yang model 
(see text) as a function of the spin polarization ( for 7 = 10. 
The exact result, obtained from the solution of the Bethe- 
Ansatz equations ((6])- ((9]), is compared with the fitting formula 
in Eq. 



Here = -1.68894, 6„ = -8.0155, = 2.74347, 
ap = -1.51457, 6/3 = 2.59864, and C/^ = 6.58046. Equa- 
tion is compared with the exact Bethe-Ansatz results 
for 7 = 10 in Fig. [TJ 

Spin-density-functional theory in the Kohn-Sham 
scheme — The ground-state spin-densities n„{z) can be 
calculated within the Kohn-Sham version of SDFT by 
solving self-consistently the effective Schrodinger equa- 
tions for single-particle orbitals 



2m dz' 



with V^l\n,]iz) = V^''\n,]iz) + V^^\n^]{z) + ^^/(z). 



<i2a,cr(z) = £a,(T Va.cr (^) (14) 



together with the closure 
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Here the sum runs over the occupied orbitals and the de- 
generacy factors ri'^'' satisfy the sum rule J2a ^^^^ ~ ^<^- 
The first term in the spin-dependent effective Kohn- 
Sham potential V-^^ is the mean-field term V-^^ = 
grD'ns{z), while the second term is the exchange- 
correlation potential defined as the functional deriva- 
tive of the exchange-correlation energy Exc[na\ eval- 
uated at the ground-state density profile, V-^c'' = 
5Ey^c[ncr\l ^na{z)\Gs- Finally, the third term is the spin- 
dependent external field Vj^^ (z) = r7ia;| z^/2 — pB<^B{z). 

Equations and provide a formally exact 

scheme to calculate n^(z), but E^c needs to be approx- 
imated. The local-spin-density approximation (LSDA) 
is known to provide an excellent description of the 
ground-state properties of a variety of inhomogeneous 
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systems [5| . In the following we employ a Bethe- Ansatz- 
based LSDA functional (BALSDA) for the exchange- 
correlation potential, 

E^c[na] ^ i?ife^^^°^[7v] (16) 

where the exchange-correlation energy e^°"^ of the homo- 
geneous Gaudin-Yang model is defined by 

exr°("-;C,ffiD) = £gs(?i,C,5id) - k("-,C) 

- £h("-, C:5id) ■ (17) 

Here 

SninX,9m) = \giDn\l-e) (18) 

is the mean-field energy. 

Illustrative numerical results for the two-body problem 
— To illustrate how this procedure based on SDFT and 
BALSDA embodies the 2kp Akp crossover we present 
some numerical results for the case N = 2, which is ex- 
actly solvable as shown in Ref. and in the Appendix of 
Ref. 0. 

For the case = 2 we choose the following form for 
the infinitesimal spin-symmetry-breaking magnetic field: 

B{z) — arctan (z/^) , (19) 

TT 



where the field strength e is reduced to zero during the 
self-consistent solution of the Kohn-Sham equations and 
^ is a parameter that we take as ^ = 0.5a||. Two impor- 
tant remarks can be made: (i) if A is small the addition 
of this infinitesimal field is irrelevant and, at the end of 
the self-consistency cycle, one reproduces the results pre- 
sented in the top and central panels of Fig. 8 in Ref. [3; 
and (ii) if A is large instead, the final result for the total 
density n{z) is very different from that resulting from the 
BALDA theory of Ref. [3. We have carefully checked that 
the final result is independent of the special form of the 
spin-symmetry-breaking field. In fact one could alterna- 
tively start the self-consistent Kohn-Sham cycle from an 
initial guess which has slightly broken spin sysmmetry 
(i.e. from n] {z) ^ n|(z)), without using any B field. 



Figure [H reports our numerical results for the total 
density profile of A^ = 2 fermions with antiparallel spins. 
The emergence of a Wigner molecule at strong coupling 
is evident and very well described by the BALSDA ap- 
proach. 
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FIG. 2; Density profile •n{z) (in units of ^) as a function of 
z/ay for two fermions with antiparallel spins at A = 1 (top) 
and A = 100 (bottom). The results of the BALSDA scheme 
proposed in this work are compared with the exact solution of 
the two-body problem and with the BALDA results of Ref. 
At A = 1 the BALSDA result coincides with the BALDA one. 



